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ABSTRACT 

We  discuss  how  simultaneous  observations  by  multiple  heliospheric  imagers  can  provide  some  important 
information  about  the  azimuthal  properties  of  Coronal  Mass  Ejections  (CMEs)  in  the  heliosphere.  We 
propose  two  simple  models  of  CME  geometry  that  can  be  used  to  derive  information  about  the  azimuthal 
deflection  and  the  azimuthal  expansion  of  CMEs  from  SECCHI/HI  observations.  We  apply  these  two 
models  to  four  CMEs  well-observed  by  both  STEREO  spacecraft  during  the  year  2008.  We  find  that  in 
three  cases,  the  joint  STEREO-A  and  B  observations  are  consistent  with  CMEs  moving  radially  outward. 
In  some  cases,  we  are  able  to  derive  the  azimuthal  cross-section  of  the  CME  fronts,  and  we  are  able  to 
measure  the  deviation  from  self-similar  evolution.  The  results  from  this  analysis  show  the  importance  of 
having  multiple  satellites  dedicated  to  space  weather  forecasting,  for  example  in  orbits  at  the  Lagrangian 
L4  and  L5  points. 

Subject  headings:  scattering  —  MHD  -  Sun:  corona  —  Sun:  coronal  mass  ejections  (CMEs) 


1.  INTRODUCTION 

With  the  launch  of  the  Solar  Mass  Ejection  Imager 
(SMEI)  and  the  Solar  Terrestrial  Relations  Observa¬ 
tory  (STEREO)  in  2003  and  2006,  respectively,  Coro¬ 
nal  Mass  Ejections  (CMEs)  can  now  be  observed  re- 
motely  all  the  way  to  the  Earth  (jHarrison  et  al.l  [20091). 
and  their  properties  compared  with  in-situ  observations  at 
1  AU  (jDavis  et  al.ll200flHRouilla.rd  et  al.ll2009lh  Although 
the  Heliospheric  Imagers  (HI-1  and  HI-2,  see  lEvles  et  al.l 
l2009t) .  part  of  the  Sun  Earth  Connection  Coronal  and  He- 
liospheric  Investigation  (SECCHI)  suite  onboard  STEREO 
(iHoward  et  al .Il2002ll.  enable  the  tracking  of  CMEs,  they 
cannot  provide  directly  the  CME  position  because  what  is 
observed  is  the  Thomson  scattered  signal  integrated  over  a 
line-of-sight.  Therefore,  the  CME  position  and  kinematics 
have  to  be  derived  from  a  variety  of  models,  which  make 


assumptions  regarding  the  CME  shape  and  its  direction 
of  propagation  (fVourlida,s  fc  Howard  |2006|:  Sheelev  et  alJ 


Il999t  IHoward  fc  Tanninl  2009!) .  In  |Lugaz_e^dJ  (20091 
using  real  and  simulated  data  from  iLugaz  et  al.l  (12009a 


we  compared  some  of  the  assumptions  typically  used:  i) 
Fixed-^,  where  it  is  assumed  that  the  His  track  a  sin¬ 
gle  plasma  element  moving  radially  outward,  ii)  Point-P, 
where  the  CME  is  treated  as  an  expanding  sphere  cen¬ 
tered  on  the  Sun,  and  iii)  a  newly  proposed  harmonic  mean 
(HM)  model,  where  the  CME  cross-section  is  a  circle  an¬ 
chored  at  the  Sun  and  its  diameter  corresponds  to  the  har¬ 
monic  mean  of  the  positions  derived  by  the  Fixed-</>  and 
Point-P  approximations.  We  found  that  the  HM  model 
gives  the  best  results  for  a  wide  CME  observed  by  one 
spacecraft  on  the  limb,  while  the  Fixed-</>  approximation, 
which  is  the  most  widely  used,  gives  comparable  results 
up  to  approximative^  0.5  AU. 

Most  of  the  work  so  far  on  SECCHI/HI  data  has  focused 


on  observations  by  a  single  spacecraft,  whereas  stereo¬ 
scopic  CME  observations  by  the  SECCHI  coronagraphs 
(COR-1  and  COR-2)  have  resulted  in  a  number  of  studies 
to  improve  our  understanding  of  CME  physical  param¬ 
eters.  For  example,  iThernisien  et  al.l  (j 2 00911  used  multi¬ 
spacecraft  observations  to  determine  the  CME  direction, 
speed  and  orientation,  and  IColaninno  &;  Vourlidas!  (|2009T) 
used  multi-spacecraft  observations  to  determine  the  CME 
direction  and  mass.  In  ILugaz  et  al.l  ('20051.  we  showed 
how  CME  images  made  by  wide-angle  white-light  im¬ 
agers  from  different  viewpoints  can  help  in  determining 
the  CME  direction  of  propagation  in  the  heliosphere.  Re¬ 
cently,  |^^oT^^owardJ2002)_pro£Osed  a  technique  sim¬ 
ilar  to  that  of  IThernisien  et  akl  ((2009)  but  which  is  ex¬ 
panded  to  the  HI  fields-of-view  and  with  a  full  treatment 
of  the  Thomson  scattering.  The  authors  use  a  geometri¬ 
cal  model,  which  is  fitted  visually  to  joint  observations  by 
the  two  STEREO  spacecraft  to  derive  the  CME  aspect, 
density  structure  and  orientation.  iLiu  et  all  (I2010T) .  in  a 
recent  study,  calculated  the  CME  central  longitude  by  tri¬ 
angulation  in  the  COR  and  HI  field-of-views  for  one  CME 
event  and  obtained  good  results  up  to  approximatively 
0.5  AU. 

In  this  article,  we  discuss  how  simultaneous  observa¬ 
tions  of  a  CME  from  the  two  STEREO  spacecraft  can 
be  used  to  derive  the  CME  central  longitude  or  its  ra¬ 
dius  of  curvature,  in  addition  to  its  radial  distance.  To 
do  so,  we  propose  two  simple  models  of  the  CME  de¬ 
scribed  in  section  [2j  Contrary  to  triangulation  techniques, 
we  do  not  assume  that  the  His  are  able  to  track  the  ex¬ 
act  same  feature  for  both  spacecraft,  but  that  they  ob¬ 
serve  different  parts  of  the  same  structure.  In  the  first 
model,  the  CME  cross-section  is  treated  as  an  expanding- 
propagating  circle  attached  to  the  Sun  as  proposed  sepa- 
rately  inlWebb  et  alJ  (l2009fl .  [Howard  fe  Tappinl  (|2009T)  and 
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Fig.  1. —  Sketches  of  the  the  two  models  used  and  described  in  this  work,  illustrating  two  ways  to  explain  asymmetric  observations.  On 
the  left,  the  CME  is  modeled  as  a  sphere  connected  to  the  Sun  with  a  varying  direction  of  propagation  </>.  On  the  right,  the  CME  is  modeled 
as  a  sphere  on  a  fixed  direction  of  propagation  (p  =  —21°  but  with  a  radius,  Ri  varying  with  time.  The  sketches  are  to  scale  except  for  the 
size  of  the  Sun,  the  Earth  and  the  spacecraft  and  correspond  approximatively  to  the  geometry  on  April  26,  2008.  The  angles  are  ea  =  39° 
and  eb  =  18°.  The  varying  cp  model  is  with  <p  =  —39°  and  R  =  135  Rq  and  the  varying  radius  model  with  Ri  =  35  Rq  and  R2  =  90  Rq. 


iLugaz  et  al.l  (|2009bf).  but  with  a  varying  central  position. 
In  the  second  model,  the  CME  has  a  fixed  direction  but 
it  is  not  attached  to  the  Sun  and  its  diameter  can  vary 
freely.  We  apply  these  two  models  to  the  analysis  of  4 
CMEs  observed  in  2008  in  section  [3]  The  conclusions  of 
this  investigation  are  drawn  in  section  [4j 

2.  DERIVATION  OF  CME  CENTRAL  LONGITUDE  AND 

RADIUS  FROM  MULTIPOINT  OBSERVATIONS 

2.1.  Model  1:  CME  Central  Longitude 

CMEs  propagating  between  the  STEREO-A  and  B 
spacecraft,  i.e.  CMEs  propagating  approximatively  to¬ 
wards  Earth,  can  be  imaged  to  large  elongation  angles  by 
the  His  onboard  both  STEREO  spacecraft.  Taking  into 
account  the  difference  in  spacecraft  heliocentric  distances 
(i dA  ~  0.95  AU,  dB  ~  1-05  AU),  the  fact  that  the  His 
observe  the  same  CME  at  the  same  time  at  different  elon¬ 
gation  angles  can  give  us  information  about  the  CME  cen¬ 
tral  position.  Such  information  cannot  be  deduced  using 
the  Point-P  approximation,  because  it  assumes  a  too  sim- 
plistic  CME  geometry.  Direct  triangulation  can  be  done 
dLiu  et  ahll2010H.  but  only  under  the  assumption  that  both 
STEREO  spacecraft  observe  the  exact  same  plasma  ele¬ 
ment,  which  is  unlikely  to  be  the  case  for  most  CMEs 
at  large  elongation  angles.  For  some  CMEs,  direct  tri¬ 
angulation  might  not  give  realistic  results.  We  suspect 
that  it  is  the  case  for  the  April  26,  2008  CME,  because  it 
appears  as  an  halo  CME  in  STEREO-B.  The  best  fit  of 
the  method  of  iRouillard  et  al.l  (120081)  for  this  CME  gives 
(j>A  =  —33.5  ±  18.0°  and  (j>B  =  2.1  ±  6.5°  for  STEREO-A 
and  B  respectively;  these  two  results  are  not  consistent 
with  each  other.  However,  one  can  use  the  expanding- 
propagating  bubble  approximation  of  iHoward  fe  Tappir] 
(200^)  and  iLueaz  et  al.l  (l2009bll  with  two  parameters  to 
analyze  the  data  from  the  two  spacecraft  simultaneously; 
for  each  pair  of  position  angles  (for  example  90  and  270), 
the  derived  parameters  are  the  diameter  of  the  circular 
front  and  its  direction  of  propagation  (see  left  panel  of 


Figure  1).  The  main  assumptions  of  this  model  are  that 
the  CME  cross-section  is  circular  and  that  the  CME  is 
anchored  at  the  Sun. 

Assuming  simultaneous  observations  by  STEREO-A 
and  STEREO-B ,  one  can  write  an  expression  for  the  di¬ 
ameter  of  the  sphere  for  each  of  the  spacecraft  following 
ILugaz  et  all  (2009bjl: 

a. = ,,  a) 

1  +  sin (/3a  ~  (p  +  eA ) 

where  <j>  is  the  direction  of  propagation  of  the  CME  away 
from  the  Sun-Earth  line  (defined  positive  from  B  to  A) 
and  Pa  is  the  angular  separation  (defined  as  a  positive 
number)  of  the  spacecraft  with  the  Sun-Earth  line.  The 
same  relation  applies  for  spacecraft  B  by  replacing  <f>  by 
—<j).  Solving  for  Ra  =  Rb  and  regrouping  the  terms  in  <f> 
gives: 

p-r 


=  arcsm 


Q 


(2) 


with 

P  =  dB  sin  eB/{dA  sin  ca), 

Q  =  y/ P ^  T  2 P  cos (fdB  T  Pa  A  eB  -T  ca)  T  1?  and 

Ps'm(pA  +  £a)  -  sui(/3b  +  eB) 

tan  a  =  — - 7— - r - 7- - 7 . 

P  cos  IP  a  +  £a)  +  cos(/3b  +  £B) 


The  diameter  of  the  sphere  can  then  be  simply  derived 
from  equation  (1). 

2.2.  Model  2:  CME  Radius 

Another  totally  independent  way  to  explain  the  differ¬ 
ent  elongation  measurements  of  the  two  spacecraft  is  to 
consider  that  it  is  due  to  the  expansion  of  the  CME  front 
(see  right  panel  of  Figure  1).  In  this  model,  the  CME  is 
not  anchored  at  the  Sun  and  its  radius  is  one  of  the  derived 
parameters.  To  compute  the  CME  radius,  we  have  to  con¬ 
sider  the  direction  of  propagation,  <j>,  as  a  known  quantity. 
It  can  be,  for  example,  obtained  from  the  COR-2  fitting 
done  by  iThernisien  et  al.l  (|2009T) .  The  main  assumptions 
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Fig.  2. —  April  26,  2008  CME  as  observed  by  SECCHI-A/ COR-2  {left)  and  HI-1  on  STEREO-A  {midde)  and  STEREO-B  {right).  The  His 
images  are  taken  33  seconds  apart,  illustrating  the  asymmetric  nature  of  the  observations.  The  yellow  dots  illustrate  the  positions  tracked  by 
the  NRL  measurements  (top  and  center  of  the  ejecta). 


of  this  model  are  the  circular  cross-section  and  the  absence 
of  heliospheric  deflection. 

Noting  R\  as  the  radius  of  the  CME,  R2  as  the  heliocen¬ 
tric  distance  of  the  center  of  the  CME,  Rc  as  the  distance 
of  the  center  of  the  CME  to  the  observer,  and  6a  —  7  as 
the  angle  between  the  Sun,  STEREO-A  and  the  center 
of  the  CME,  the  law  of  cosines  and  the  law  of  sines  give, 
respectively: 

Rc  =  d2A+  R2  —  2(1aR2  cos(/3a  —  <t>) 

R2  sm(/3A  -<j>)=Rc  sin(eA  -  7) 

But  Rc  sin  7  =  Ri  and  Rc  cos  7  =  1 J R2;  —  Rf  and  the 
second  equation  can  be  rewritten  as: 

R2  sin(/3J4  —  <j>)  =  y/ R^,  —  R\  sin  ca  —  R\  cosca 

which  can  be  combined  with  the  first  equation  into  the 
following  simple  equation: 

R2  sin(/?A  ~4>)  +  R\  cos  £a 

=  sin  €a \J  d?A  +  R2  -  R\  -  2cIaR2  cos  {Pa  -  4>),  (3) 

with  the  equation  for  spacecraft  B  obtained  by  replac¬ 
ing  <j>  by  —cj).  By  solving  equations  (3)  for  STEREO- 
A  and  STEREO-B  simultaneously,  R\  and  R2  can  be 
uniquely  derived.  The  CME  half-angle  is  given  by  0  = 
arctan(i?i/i?2)-  It  should  be  noted  that  for  CMEs  prop¬ 
agating  directly  towards  Earth  (</>  ~  0),  this  model  al¬ 
ways  predicts  that  the  HI  instruments  onboard  the  two 
STEREO  spacecraft  will  observe  the  exact  same  time- 
elongation  profile,  which  is  a  direct  consequence  of  the 
assumption  of  a  circular  front.  Also,  if  ca  <  £b  at  one 
time  but  cb  <  ca  at  another  time,  there  is  no  solution 
with  I?i  >  0  and  R2  >  0,  and  the  model  cannot  be  ap¬ 
plied.  In  these  cases,  the  CME  deformation,  deflection 
and  the  deviation  from  circular  of  its  cross-section  must 
be  taken  into  account. 

The  main  hypothesis  of  these  two  models  is  that  the  HI 
instruments  onboard  STEREO-A  and  B  do  not  observe 
the  same  plasma  element,  as  is  likely  to  be  the  case  in  the 
COR  fields-of-view  and  as  is  assumed  in  iLiu  et  all  (|2010l). 
The  assumption  here  is  that  the  CME  front  is  locally  circu¬ 
lar  and  that  it  is  “projected”  onto  the  HI  fields-of-view  at 


the  elongation  angle  corresponding  to  the  tangent  to  this 
front.  Similar  assumptions  have  been  made  bvlWood  et  al.l 
(12009!)  and  iHoward  fc  Tanninl  (|2009l)  to  derive  CME  po¬ 
sitions  from  observations  at  large  elongation  angles  from 
one  spacecraft. 

3.  DATA  ANALYSIS 

3.1.  Data  Selection 

For  these  new  methods  to  be  applied,  the  studied  CMEs 
have  to  propagate  between  the  two  STEREO  spacecraft, 
and  they  may  appear  as  halo,  partial  halo  or  wide  CMEs 
in  LASCO  field-of-view.  To  find  potential  CMEs  to  study, 
we  look  at  all  instances  of  CMEs  with  an  apparent  width 
larger  than  100°  in  LASCO  field-of-view  from  January 
2008  when  the  spacecraft  separation  reached  45°  until  May 
2009.  We  found  17  CMEs,  5  of  which  were  observed  si¬ 
multaneously  by  the  HI  instruments  onboard  STEREO-A 
and  B ,  they  are  the  April  26,  June  2,  August  30,  Novem¬ 
ber  3  and  December  12,  2008  CMEs.  The  November  3  is 
associated  with  an  instance  of  CME-CME  interaction  and 
we  exclude  it  from  the  current  study  (jKilpua  et  al.ll2009h . 
Although  the  August  30  CME  also  interacts  with  a  pre¬ 
ceding  ejection,  the  preceding  ejection  appears  to  be  small 
enough  in  size  so  that  it  did  not  influence  too  strongly 
the  propagation  of  the  overtaking  CME  and,  therefore,  we 
keep  this  CME  in  this  study.  The  December  12  CME  had 
two  bright  fronts  observed  in  both  spacecraft  and  we  ana¬ 
lyzed  3  different  datasets  for  the  April  26  CME,  resulting 
in  a  total  of  7  analyzed  pairs  of  datasets.  We  first  ap¬ 
ply  the  two  models  to  a  detailed  case  study  for  the  April 
26,  2008  CME,  before  presenting  the  results  for  the  other 
CMEs. 

3.2.  April  26  CME 

On  April  26,  2008  at  1425UT,  SECCHI-A/COR-1  ob- 
served  an  eastern  limb  CME.  This  event  was  studied 
by  Thernisien  et  al.l  (|2009f)  and  IColaninno  fc  Vourlidasi 
d2009|)  in  the  COR-1  and  COR-2  fields-of-view  and  by 
IWood  fc  Howardl  (|2009l)  in  the  entire  SECCHI  field-of- 
view.  Images  of  the  CME  in  COR2-A,  HI1-A  and  HI1-B 
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Fig.  3. —  Top  Left:  Central  longitude,  <f>,  as  derived  from  the  varying-^  model  (model  1).  Top  Right:  CME  radius,  R\  (solid),  and  half-angle, 
0  (dash-dot),  as  derived  from  the  varying-radius  model  (model  2).  Bottom  Left:  Elongation  measurements  of  the  CME,  showing  the  RAL 
measurement  and  the  NRL  measurements  for  the  top  and  the  center  of  the  CME.  Bottom  Right:  Error  for  the  RAL  measurements  assuming 
the  elongation  angles  are  measured  with  a  precision  of  15  pixels. 


are  shown  in  Figure  2.  To  do  our  analysis,  we  used  elon¬ 
gation  measurements  for  PA  90  and  270  obtained  from 
the  Naval  Research  Laboratory  (hereafter,  NRL  measure¬ 
ments)  from  J-maps  including  COR  and  HI  observations 
and  from  the  Rutherford  Appleton  Laboratory  (hereafter, 
RAL  measurements)  from  J-maps  derived  from  Hl-only 
observations  along  the  ecliptic.  The  time-elongation  plots 
are  shown  on  the  bottom  left  panel  of  Figure  3.  The 
tracked  front  corresponds  to  the  top  of  the  ejecta  for  the 
RAL  and  the  first  of  the  NRL  measurements  and  to  the 
center  of  the  ejecta  for  a  second  set  of  NRL  measurements. 
Two  difficulties  with  applying  the  models  described  above 
to  real  data  are  that  the  elongation  measurements  must 
be  derived  accurately  and  that  the  boundary  of  the  same 
structure  must  be  tracked  in  STEREO-A  as  well  as  in 
STEREO-B.  This  last  condition  may  not  always  be  ful¬ 
filled  in  HI-2  where  the  cause  of  bright  fronts  is  sometimes 
hard  to  establish  (e.g.,  seelLugaz  et  al.l[2008l). 

We  analyzed  the  time-elongation  data  for  STEREO-A 
and  B  for  this  CME  with  the  two  methods  proposed  above. 
Using  data  from  two  different  groups  allow  us  to  quantify 
the  errors  associated  with  the  manual  determination  of  the 
elongation  angles.  The  results  of  our  analysis  are  shown 
in  Figure  3.  In  the  bottom  right  panel  of  this  Figure,  we 
show  the  error  bars  for  model  1  (black)  and  model  2  (red) 
assuming  the  elongation  measurements  are  made  with  a 
precision  of  15  pixels  corresponding  to  an  uncertainty  of 
±.15°  and  ±.5°  in  HI-1  and  HI-2  helds-of-view,  respec¬ 
tively.  The  error  in  direction  (model  1)  is  typically  ±2°  for 
both  instruments  and  the  error  in  the  CME  radius  (model 


2)  is  ±1  Rq  in  HI-1  and  increases  to  ±2.5  Rq  when  the 
CME  is  observed  in  both  HI-2  simultaneously. 

According  to  model  1,  the  CME  is  deflected  towards 
the  East  (i.e.  away  from  Earth)  reaching  a  near-constant 
central  longitude  of  -35±2°  at  about  50  Rq,  the  CME 
continues  to  be  deflected  towards  the  East  with  a  rate  of 
about  -3.5°/day  from  this  distance  on  (see  top  left  panel 
of  Figure  3).  The  initial  angle  is  ~  —20°  comparable 
to  the  angle  of  -21°  derived  by  iThernisien  et  alJ  (|2009t) 
based  on  COR2  data.  The  central  position  of  the  CME 
using  only  STEREO-A  data  and  the  Fixed-/)  procedure 
of  iRouillard  et  ah  ( 200811  can  be  estimated  at  -33.5±18°, 
while  Wood  fe  Howardl  ( 20091)  reports  a  best-fit  value  of 
-26.5°.  Our  value  of  -35  ±  2°  appears  in  relatively  good 
agreement  with  these  values. 

For  model  2,  we  assume  the  fixed  direction  of  propa¬ 
gation  />  =  —21°  from  IThernisien  et  al.l  (|2009l).  Then, 
the  CME  radius  monotonically  increases  from  a  value  of 
8±3  Rq  at  20  Rq  to  26±1  Rq  at  80  Rq,  and  to  about 
37  Rq  at  140  Rq.  The  corresponding  CME  half-angle  is 
equal  to  25  ±2°  from  60  Rq  to  120  Rq  with  a  general  rate 
of  shrinking  of  4° /day  until  the  end  of  the  measurements 
(see  top  right  panel  of  Figure  3). 

In  the  first  40  Rq,  the  results  using  the  different  datasets 
are  inconsistent  with  each  other  (see  top  panels  of  Fig¬ 
ure  3).  This  is  because  the  CME  appears  in  STEREO-B 
as  an  halo  and  in  STEREO-A  as  a  limb  CME,  and  we 
are  not  able  to  identify  the  same  front  in  the  observations 
from  the  two  spacecraft.  In  fact,  the  CME  is  first  visi¬ 
ble  in  HI-B  (at  ~  4°)  when  it  reaches  a  distance  of  about 
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Fig.  4. —  Left:  Central  longitude,  </>,  as  derived  with  the  varying-^  model  (model  1)  for  the  June  2  and  August  30  CMEs  and  for  the  two 
bright  features  associated  with  the  December  12  CME.  Right :  CME  radius,  R\,  as  calculated  for  model  2  for  the  June  2  CME  and  the  second 
front  of  the  December  12  CME. 


32  Rq  from  the  Sun,  corresponding  to  observations  around 
8°  elongation  in  HI- A  (see  bottom  left  panel  of  Figure  3). 
Before  this  time,  since  there  is  no  measurement  of  the 
CME  in  HI-B,  we  use  the  elongation  angles  measured  by 
STEREO-B  in  COR-2  for  the  NRL  data.  We  find  it  im¬ 
possible  to  separate  the  top  of  the  ejecta  from  the  shock 
front  and  piled-up  mass  in  COR2-B  and  early  on  in  HI1- 
B.  Therefore,  we  are  not  able  to  identify  the  same  front 
between  STEREO-A  and  STEREO-B  when  the  CME  is 
within  40  Rq  from  the  Sun,  and  the  results  are  inconsis¬ 
tent  between  the  RAL  and  NRL  measurements.  For  the 
center  of  the  ejecta,  we  are  able  to  identify  the  common 
feature  in  both  STEREO-A  and  B  and  the  models,  in  this 
case,  give  a  more  realistic  evolution  of  <f>  and  R\.  This 
is  because  in  STEREO-B ,  where  the  CME  appears  as  an 
halo,  it  is  still  relatively  easy  to  separate  the  center  of  the 
ejecta  from  the  other  density  features.  It  should  be  noted 
that  the  center  of  the  ejecta  might  be  of  relatively  small 
angular  extent  in  the  azimuthal  direction,  in  which  case, 
using  direct  triangulation  as  in  iLiu  et  all  (1201081  is  more 
adapted. 

We  compared  the  CME  distance  as  derived  with  these 
two  models,  with  the  HM  model  proposed  in  iLugaz  et  al.l 
(|2009bH  using  only  STEREO-A  data.  All  methods  yield 
the  same  distances  for  the  CME  within  10%.  However, 
from  the  two  new  models,  it  is  possible  to  predict  a 
hit/miss  at  different  spacecraft’s  positions  in  the  helio¬ 
sphere  based  on  the  azimuthal  properties  of  the  CME. 
Based  on  the  height-time  profile  in  the  heliosphere  us¬ 
ing  the  second  model,  we  fit  the  data  to  a  CME  with  a 
constant  speed  of  534  km  s-1  -to  compare  to  the  final 
speed  of  543  km  s-1  derived  bvlWood  fe  Howardl  (12009H- 
.  Additionally,  we  calculate  the  CME  half-angle,  9 ,  with 
the  fixed  rate  of  — 4°/day,  using  the  best-fit  formula  of 
9  =  25°  —  4°  x  (f  —  to))  where  t  is  the  time  in  day  since  R 
=  April  27  1200UT.  With  these  parameters,  the  model  cor¬ 
rectly  predict  that  the  CME  does  not  hit  ACE  but  it  pre¬ 
dict  that  the  CME  hit  STEREO-B  at  0300UT  on  04/30. 
In-situ  observations  by  ACE  and  the  two  STEREO  space¬ 
craft  show  that  only  STEREO-B  detected  the  passage  of 
an  iCME  from  1530UT  on  04/29  to  0700UT  on  04/30, 
which  translates  into  an  error  of  about  11  hours  for  the 


arrival  time  for  our  model. 


3.3.  June  2,  August  30  and  December  12  CMEs 

The  June  2,  2008  CME  was  included  in  the  study  by 
iThernisien  et  al.l  (200911  and  studied  in  iRobbrecht  et  alJ 
(|2009l).  while  the  December  12,  2008  CME  has  been  ana¬ 
lyzed  bv  I  Davis  et  all  (|2009h  and  ILiu  et  al~l  (|2010H.  For  our 
study,  we  used  the  data  available  from  the  Rutherford  Ap¬ 
pleton  Laboratory  website.  In  the  left  panel  of  Figure  4,  we 
show  the  variation  of  the  central  longitude  of  the  CMEs  fol¬ 
lowing  model  1  for  the  4  datasets.  In  Table  1,  we  compare 
our  results  with  information  available  from  the  flare  loca¬ 
tion,  with  the  direction  of  propagation  as  calculated  with 
the  procedure  of  iRouillard  et  all  (l2008il  and  with  other 
published  analyses.  It  should  be  noted  that  the  method 
of  ILiu  et  alJ  (20ld)  assumes  that  Ja  =  ds,  whereas  we 
use  the  real  values  of  the  spacecraft  heliocentric  distances 
(which  we  assume  constant  over  the  duration  of  an  event). 
In  our  experience,  doing  so  is  required  to  obtain  consistent 
values,  especially  at  large  elongation  angles.  For  example, 
analyzing  the  last  pair  of  elongation  angles  from  our  mea¬ 
surements  for  the  second  front  of  the  December  12  event, 
the  method  of  iLiu  et~.'.  (12010)  gives  —5°  with  Ja  =  ds, 
but  using  the  exact  values  of  Ja  and  dg  shifts  the  derived 
longitude  from  -5°  to  13°.  This  corrected  value  is  closer  to 
the  last  value  of  27°  obtained  with  model  1  than  the  value 
of  —5°  reported  bvlLiu  et  al.l  (1201011. 

We  find  that  the  June  2  CME  and  the  two  fronts  from 
the  December  12  CME  propagate  close  to  radially  out¬ 
ward  until  about  140  Rq  (0.65  AU),  while  the  August  30 
CME  shows  a  deflection  towards  the  East  of  about  30° 
during  its  propagation  and  crosses  the  Sun-Earth  line.  As 
noted  in  section  2.2,  the  second  model,  because  it  assumes 
a  fixed  direction  of  propagation  cannot  be  used  for  CMEs 
propagating  exactly  towards  Earth  (</  ~  0°)  nor  for  CMEs 
which  cross  the  Sun-Earth  line  during  their  propagation, 
therefore  we  only  applied  this  model  with  the  measure¬ 
ments  from  the  June  2  CME  as  well  as  the  second  front  of 
the  December  12  CME,  which  appear  to  propagate  away 
from  the  Sun-Earth  line.  We  used  directions  of  propaga¬ 
tion  of  -15°and  20°  for  these  features,  respectively.  From 
model  1,  the  June  2  CME  appears  to  move  on  a  quasi- 
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CME 

Method  &  Instruments  (Reference) 

Phi  Estimate 

Apr.  26 

Flare 

-9U 

Apr.  26 

Visual  COR2s  (1) 

-21° 

Apr.  26 

Visual  STEREOs  (2) 

-28° 

Apr.  26 

Mass  COR2s  (3) 

-48° 

Apr.  26 

Fixed-  (f>  HI- A 

-33.5  ±  18° 

Apr.  26 

Fixed-  4>  HI-B 

2.1  ±  6.5° 

Apr.  26 

Model  1 

-17°  (at  20  Rq)  to  -37.5°  (at  130  Rq) 

June  2 

Visual  COR2s  (1) 

-37° 

June  2 

Fixed-  <fi  HI- A 

-24.2  ±  5° 

June  2 

Fixed- ^  HI-B 

20.9  ±  11.5° 

June  2 

Model  1 

-17  ±2.7° 

Aug.  30 

Fixed-0  HI-A 

-16  ±  17° 

Aug.  30 

Fixed- (f>  HI-B 

19.2  ±  10.5° 

Aug.  30 

Model  1 

10  °  (at  20  Rq)  to  -20°  (at  140  Rq) 

Dec.  12  Front  1 

Fixed- 0  HI-A 

-14.7  ±  13° 

Dec.  12  Front  1 

Fixed- ^  HI-B 

12.6  ±  6.5° 

Dec.  12  Front  1 

Triangulation  (4) 

0  ±  5° 

Dec.  12  Front  1 

Model  1 

10  ±  10° 

Dec.  12  Front  2 

Fixed-^  HI-A 

8.3  ±  4.5° 

Dec.  12  Front  2 

Fixed- 4>  HI-B 

-1.5  ±  7° 

Dec.  12  Front  2 

Triangulation  (4) 

5  ±  3°  (up  to  70  Rq)  -3  ±5°  (after) 

Dec.  12  Front  2 

Model  1 

20  ±  7° 

Table  1 

Direction  of  propagation  for  the  four  CMEs  from  this  work  compared  to  other  methods.  References:  1 
Thernisien  et  al.  (2009),  2  Wood  &  Howard  (2009),  3  Colaninno  &  Vourlidas  (2009),  4  Liu  et  al.  (2010). 


radial  trajectory  at  about  —15°,  while  the  second  front 
of  the  December  12  CME  appears  to  propagate  on  a  di¬ 
rection  of  about  20°  with  respect  of  the  Sun-Earth  line; 
hence,  we  choose  these  numbers  for  the  fixed  directions  of 
propagation.  The  results  of  the  second  model  are  shown 
in  the  right  panel  of  Figure  4. 

The  June  2  observations  can  be  explained,  using  model 
2,  by  a  CME  which  propagates  on  a  fixed  radial  trajectory 
and  whose  radius  in  the  ecliptic  increases  more  slowly  than 
self-similarly,  the  CME  half-angle  decreasing  from  45°  to 
35°  (see  Figure  4).  Alternatively,  it  can  be  explained  by 
model  1  as  a  CME  whose  direction  of  propagation  is  de¬ 
flected  by  about  5°  towards  the  East  in  about  0.45  AU. 
Both  these  explanations  appear  physically  realistic,  involv¬ 
ing  a  propagation  close  to  radially  outward  and  an  evolu¬ 
tion  close  to  self-similar,  and  it  is  likely  that  the  evolution 
of  the  June  2  CME  is  a  combination  of  these  two  results, 
with  a  limited  eastward  deflection  and  a  small  shrinking 
of  the  CME  front.  The  observations  from  the  August  30 
CME  cannot  be  explained  by  model  2,  because  the  CME 
appears  first  farthest  in  STEREO-B  before  appearing  far¬ 
thest  in  A  (as  noted  in  2.2,  this  causes  model  2  to  be 
inapplicable).  Using  model  1,  it  corresponds  to  a  CME 
being  deflected  by  30°  towards  the  East  in  0.5  AU.  The  2 
features  from  the  December  12  CME  appear  to  propagate 
close  to  the  Sun-Earth  line  on  near-radial  trajectories  but 
they  cannot  be  simply  analyzed  with  model  2,  either  be¬ 
cause  the  fronts  propagate  too  close  from  the  Sun-Earth 
line  or  because  the  model’s  assumption  of  observing  the 
tangents  to  a  circular  CME  cross-section  is  not  correct  for 
this  CME.  It  appears  from  both  our  study  and  that  by 
iLiu  et  abl  (|2010fl  that  the  first  front  of  this  CME  propa¬ 


gates  about  5-10°  East  of  the  second  front. 

4.  DISCUSSION  AND  CONCLUSIONS 

In  this  article,  we  propose  two  models  to  derive  informa¬ 
tion  about  the  azimuthal  properties  of  CMEs  from  multi¬ 
spacecraft  observations  in  the  heliosphere  using  simple  ge¬ 
ometrical  considerations.  These  models  can  be  used  to 
derive  the  CME  radius  or  central  position  from  SECCHI 
observations  without  having  to  rely  on  any  extra  human 
judgement  or  fitting  procedure,  after  obtaining  the  time- 
elongation  measurements.  The  main  hypothesis  of  these 
two  models  is  that  the  HI  instruments  onboard  STEREO- 
A  and  B  do  not  observe  the  same  plasma  element,  as  is 
likely  to  be  the  case  in  the  COR  fields-of-view  and  as  is 
assumed  inlLiu  et~aLl  (|2010|) .  The  assumption  here  is  that 
the  CME  front  is  locally  circular  and  that  it  is  viewed  in 
the  His  at  the  elongation  angle  corresponding  to  the  tan¬ 
gent  to  this  front.  We  applied  these  two  new  models  to 
six  features  belonging  to  four  CMEs  observed  in  2008  by 
both  STEREO  spacecraft. 

For  five  of  the  six  studied  features,  the  April  26  (center 
and  front),  June  2  and  December  12,  2008  (front  1  and 
front  2)  CMEs,  we  find  that  the  measurements  can  be  ex¬ 
plained  as  being  from  CMEs  propagating  close  to  radially 
outward.  However,  for  two  of  the  CMEs  (the  June  2  CME 
and  both  features  from  the  December  12  CME),  we  do 
not  assume  a  radial  propagation  and  we  derive  the  CME 
central  longitude  directly  with  model  1.  We  find  that  the 
CME  central  longitude  remains  within  15°  of  radial,  with 
the  June  2  CME,  in  particular,  being  deflected  monoton- 
ically  towards  the  east  by  about  5°  in  0.5  AU.  For  these 
CMEs,  model  2,  which  assumes  radial  propagation  but 
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does  not  assume  a  self-similar  expansion,  does  not  provide 
additional  information.  Also,  the  two  features  associated 
with  the  December  12  CME  are  found  to  be  separated  by 
about  10°  at  all  time,  which  was  also  found  using  direct 
triangulation  bvlLiu  et~aT.'  tl20ldl. 

For  the  April  26,  2008  CME,  the  model  assuming  self¬ 
similar  expansion  (model  1)  appears  to  fail,  because  it  pre¬ 
dicts  a  large  deflection  of  the  CME  which  is  not  phys¬ 
ically  expected  nor  confirmed  by  in-situ  measurements. 
Self-similar  expansion  of  CMEs  in  the  heliosphere  has 
long  been  assumed  (IXue  et  al.ll20f)5t  iKrall  fc  St.  CvifeOQfi: 

IWood  fe  Howard!  I2009H .  but  in  general,  it  has  not  been 
tested  with  heliospheric  observations.  We  analyze  this 
CME  with  model  2  ,  where  self-similar  expansion  is  not 
assumed  and  using  the  direction  of  propagation  derived 
by  iThernisien  et  alJ  (2009)  ■  We  derive  the  change  over 
time  of  the  CME  cross-section,  which  is  found  to  decrease 
with  a  rate  of  about  4°/day.  The  cause  of  this  decreasing 
cross-section  has  to  be  further  studied,  but  it  shows  that 
the  CME  expands  slower  than  self-similarly.  This  may,  for 
example,  reflect  a  variation  of  the  CME  radius  of  curvature 
in  the  ecliptic  plane.  This  result  is  found  for  solar  mini¬ 
mum  conditions  when  the  background  solar  wind  is  more 
simply  structured.  For  the  other  CME  (August  30,  2008) 
the  measurements  can  only  be  explained,  with  model  1, 
by  a  a  deflection  towards  the  East  of  30°  in  0.5  AU  of  the 
CME.  It  is  also  possible  that  an  instance  of  CME- CME 
interaction  locally  deformed  this  CME  front  away  from 
circular  which  would  render  our  models  less  accurate. 

In  this  study,  we  ignore  the  effect  of  the  Thomson  sphere 
and  the  angular  dependence  of  the  Thomson  scattering  in¬ 
tensity.  To  apply  the  new  analysis  techniques  presented 
here,  the  CME  has  to  propagate  between  the  STEREO- 
A  and  STEREO-B  spacecraft.  For  this  study  focusing 
on  CMEs  observed  in  2008,  this  means  that  the  CMEs 
propagate  within  40°  from  the  Sun-Earth  line,  and  we  be¬ 
lieve  that  ignoring  the  difference  in  Thomson  scattering 
between  the  two  spacecraft  is  justified  as  a  first  approxi¬ 
mation.  However,  as  the  two  STEREO  spacecraft  continue 
to  separate,  CMEs  propagating  farther  away  from  the  Sun- 
Earth  line  will  be  observed  simultaneously  by  both  space¬ 
craft  and  a  more  thorough  analysis  of  the  effect  of  the 
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two  different  Thomson  spheres  should  be  made.  Addition¬ 
ally,  effects  associated  with  the  interaction  of  a  CME  with 
the  Thomson  sphere  have  been  found  to  happen  at  elonga¬ 
tion  angles  greater  than  45°  (e.g.  iManchester  et  al.l  (l2068h. 
iLugaz  et  al.l  (|2008fP.  which  is  farther  than  that  analyzed 
here. 

We  believe  that,  by  providing  a  value  of  the  CME  ra¬ 
dius  in  the  ecliptic  plane  or  the  time-dependent  variation 
of  the  central  longitude  which  can  be  simply  computed 
from  time-elongation  measurements,  these  models  could 
improve  space  weather  forecasting  of  CMEs.  However,  to 
use  such  techniques  on  a  real-world  basis,  it  might  be  use¬ 
ful  to  have  dedicated  spacecraft  making  heliospheric  obser¬ 
vations  at  fixed  locations,  for  example  from  the  Lagrangian 
points  L4  and  L5.  In  general,  we  believe  our  model  1  can 
be  used  to  derive  the  CME  central  longitude  and  its  tem¬ 
poral  variation  but  is  not  applicable  for  these  CMEs  whose 
expansion  is  not  self-similar.  When  this  is  the  case,  for 
example  because  of  interaction  with  the  structured  solar 
wind  or  with  previous  CMEs,  we  propose  a  second  model, 
which  can  quantify  the  deviation  from  self-similar,  assum¬ 
ing  radial  propagation.  Future  work  should  focus  on  the 
radius  and  central  longitude  of  a  CME  at  different  position 
angles  and  on  further  testing  and  validation  of  the  models. 
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